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Abstract: Contour ridge systems may lead to seepage that could result in serious soil erosion. Modeling soil 
erosion under seepage conditions in a contour ridge system has been overlooked in most current soil erosion 
models. To address the importance of seepage in soil erosion modeling, a total of 23 treatments with 3 
factors, row grade, field slope and ridge height, in 5 gradients were arranged in an orthogonal rotatable 
central composite design. The second-order polynomial regression model for predicting the sediment yield 
was improved by using the measured or predicted seepage discharge as an input factor, which increased the 
coefficient of determination (R?) from 0.743 to 0.915 or 0.893. The improved regression models combined 
with the measured seepage discharge had a lower P (0.007) compared to those combined with the predicted 
seepage discharge (P=0.016). With the measured seepage discharge incorporated, some significant 
(P<0.050) effects and interactions of influential factors on sediment yield were detected, including the row 
grade and its interactions with the field slope, ridge height and seepage discharge, the quadratic terms of 
the field slope and its interactions with the row grade and seepage discharge. In the regression model with 
the predicted seepage discharge as an influencing factor, only the interaction between row grade and seepage 
discharge significantly affected the sediment yield. The regression model incorporated with predicted 
seepage discharge may be expressed simply and can be used effectively when measured seepage discharge 
data are not available. 
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1 Introduction 


Land degradation caused by soil erosion is a global problem for sustainable agriculture and has 
attracted great attention (Changere et al., 1995; Ni et al., 2017; Nishigaki et al., 2017). Contour 
ridge systems, unlike flat surface and up and down ridge tillage on slopes, have the characteristic 
of water storage in the furrow between the two adjacent ridges, which leads to increased water 
infiltration and less soil erosion (Quinton and Catt, 2004; Shi et al., 2004; Brunner et al., 2008; 
Gebreegziabher et al., 2009; Grum et al., 2017). Soil erosion in contour ridge systems is more 
carefully considered in the Revised Universal Soil Loss Equation, Version 2 (RUSLE2) (USDA- 
ARS, 2008a; Liu et al., 2014a) than in other erosion models, e.g., Water Erosion Prediction Project 
(WEPP) (Flanagan and Livingston, 1995), Limburg Soil Erosion Model (LISEM) (Hessel et al., 
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2003) and Chinese Soil Loss Equation (CSLE) (Shi et al., 2013). The soil conservation benefit for 
contour ridges in RUSLE2 is presented as a subfactor to support practices in which the ridge height 
and row grade factors are quantified (USDA-ARS, 2008a). In the WEPP model, the factors of ridge 
height and row grade along with contour row spacing and contour row length are used to calculate 
the water storage in the furrow and when the overflow occurs (Flanagan and Livingston, 1995). In 
the LISEM and CSLE models, the conservation ability of the contour ridge is assigned as a constant 
in the soil loss estimations (Hessel et al., 2003; Shi et al., 2013). With the assumption of no 
infiltration and no soil erosion occurring, Razafison et al. (2012) established a model to simulate 
shallow water flow over the ridges; however, this model could not be used for soil erosion 
estimation. 

In RUSLE2, the effects of ridge geometry factors, ridge height and row grade, were carefully 
managed. The ridge height that negatively affects soil erosion determines the contouring 
effectiveness and decays as a function of the precipitation amount and interrill erosion. The relative 
row grade, which is measured as the ridge-furrow orientation to the overland flow path, was 
grouped into five classes (5%, 10%, 25%, 50% and 100%) to predict soil erosion; however, the 
prediction accuracy of additional classes is not warranted in RUSLE2. Absolute row grade, defined 
as the decrease in elevation over distance along the furrows (rise/run), gives more credit for 
contouring on various slopes compared to relative row grade (USDA-ARS, 2008a). The influence 
of field slope on soil erosion is a subject of interest in RUSLE2 and described as a concave curve 
interacting with ridge height (Liu et al., 2014a). 

In a contour ridge system, when the amount of stored rainwater exceeds the storage capacity of 
the furrow, overflow will occur, resulting in contouring failure (Griffith et al., 1990; Cui et al., 2007; 
USDA-ARS, 2008a). Being a rill-interrill erosion model, RUSLE2 did not consider soil erosion 
induced by contouring failure and the subsequent ephemeral gully erosion (USDA-ARS, 2008a). 
To address this issue, Liu et al. (2014b) used a new type of experimental box in which row grade 
and field slope could be adjusted simultaneously to qualify the effect and interaction of row grade, 
field slope, ridge height, ridge width and rainfall intensity factors on soil erosion induced by ridge 
collapse. Liu et al. (2014a) also analyzed the effects on interrill and rill erosion on row sideslopes; 
however, only two factor gradients were designed in the study and not even a simple estimated 
equation for soil erosion could be established for practical use. 

In the ridge and furrow systems, some depressions are formed when the system is constructed 
because of the microtopography (Liu et al., 2016). Concentrated rainwater with higher water levels 
(i.e., higher hydraulic gradients) in these depressions could increase water infiltration and make the 
soil saturated. When the rainwater penetrates the ridge soil, seepage may form on the row sideslope, 
especially at the foot of the ridge where the effective stress of the surface soil may be reduced, 
potentially leading to additional soil erosion (Huang and Laften, 1996; Chu-Agor et al., 2008). The 
aggravated effect of seepage on soil erosion has been the subject of recent interest (Nachshon, 2016; 
Regmi et al., 2017). Compared with the drainage regime, Nouwakpo et al. (2010) found that the 
average erodibility under a seepage regime was 5.64 times greater than that under the free drainage 
regime. Soil erosion under seepage conditions was shown to be 6 times higher for a surface under 
a hydraulic gradient of 20 cm compared with a surface drained for 7 d (Huang and Laften, 1996). 
A similar multiple was observed by Zheng et al. (2000) under run-on and runoff feed conditions. 
Nouwakpo and Huang (2012) found soil erosion was 2.1 times higher in seepage conditions. In 
addition, seepage may promote rill formation (Huang and Laften, 1996; Valentin et al., 2005; 
Nouwakpo and Huang, 2012). Nouwakpo and Huang (2012) observed that channel erosion rates 
doubled under seepage conditions. Seepage could not only increase soil erosion erodibility, but also 
transport fine particle resulting sapping, slumping or even piping (Wilson et al., 2007; Fox and 
Wilson, 2010). As the converse process to seepage, seepage weathering, defined as the weathering 
process facilitated by seepage, could also reduce the parent material's cohesion and result piping, 
suffusion and a high erosivity (Sato and Kuwano, 2015; Nachshon, 2016). The impact of seepage 
in slope failure was focused (Regmi et al., 2017) and the erosion under seepage conditions in 
contour ridging systems was investigated (Liu et al., 2015), but the impact of seepage has not been 
considered in soil erosion models. 
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Before contour failure and under drainage conditions, the influences of row grade, field slope 
and ridge height on soil erosion in contour ridge systems have been carefully considered in previous 
researches (Flanagan and Livingston, 1995; USDA-ARS, 2008a). In practice, contour ridge 
collapse on sloped land is a common phenomenon in North China (Fig. 1). In our previous studies 
(Liu et al., 2015), the influences of ridge height, row grade and field slope on runoff and soil loss 
under seepage conditions were interpreted. Here, we assumed that incorporating seepage data to 
estimate sediment yield would improve the accuracy of soil erosion model. Therefore, based on the 
dataset from Liu et al. (2015, 2016), the specific objectives of this study were to (1) model seepage 
and sediment yield with the row grade, field slope and ridge height as influencing factors; and (2) 
incorporate a seepage model with a sediment yield model to improve prediction accuracy. 


Fig. 1 Contour failure on sloped land in North China 


2 Materials and methods 


2.1 Experimental design 


According to previous studies and field investigations, 3 key factors (row grade, field slope and 
ridge height) in 5 gradients for 23 treatments were arranged in an orthogonal rotatable central 
composite design (Tables 1 and 2). This experimental design has been widely used to detect the 
effect and interaction of the key factors on dependent variables by building a second-order 
polynomial regression equation (Dominguez et al., 2010; Tang, 2010; Hadjmohammadi and Sharifi, 
2012). Relative to the full factorial design, the orthogonal rotatable central composite design allows 
a substantial reduction of the treatment number while enabling direct access to the effects of the 
items in the regression model (St-Pierre and Weiss, 2009). Based on the design method, five code 
values (—1.68, —1.00, 0.00, 1.00 and 1.68) were determined for each of these 3 factors as shown in 
Table 1 (Ding, 1986). In Table 2, treatments 1-8 were arranged orthogonally with an inherent 
repetition of these 3 factors; treatments 15-23 were the repetition treatments on the central point 
with factors at the zero code level. However, there were no repetitions in treatments 9—14; therefore, 
these treatments were replicated twice in this study. The performance of the regression model at 
the treatments of 1-14 could assess the variance explanation and at the treatments of 15-23 could 
detect the experimental errors. 


Table 1 Code gradient and values of influential factors 


Factor 
Code gradient R - 
Row grade (°) Field slope (°) Ridge height (cm) 
1.68 10.0 15.0 16.0 
1.00 8.4 13.0 14.4 
0.00 6.0 10.0 12.0 
—1.00 3.6 7.0 9.6 


—1.68 2.0 5.0 8.0 
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Table 2 Orthogonal rotatable central composite design and measured values of seepage discharge and sediment 
yield 


Factor values Experiment results 
Treatment 
number Row grade Field slope Ridge height Seepage discharge Sediment yield 

(°) (°) (cm) (L/min) (kg) 

1 8.4 13.0 14.4 0.98 1.19 
2 8.4 13.0 9.6 0.69 0.61 
3 8.4 7.0 14.4 0.87 0.33 
4 8.4 7.0 9.6 0.50 0.40 
5 3.6 13.0 14.4 0.99 2.09 
6 3.6 13.0 9.6 0.70 0.14 
7 3.6 7.0 14.4 0.97 0.14 
8 3.6 7.0 9.6 0.52 0.09 
9 10.0 10.0 12.0 0.80 3.83 
10 2.0 10.0 12.0 1.21 0.17 
11 6.0 15.0 12.0 0.54 2.09 
12 6.0 5.0 12.0 0.74 0.17 
13 6.0 10.0 16.0 1.26 0.77 
14 6.0 10.0 8.0 0.43 2.03 
15 6.0 10.0 12.0 0.51 2.20 
16 6.0 10.0 12.0 0.63 3.48 
17 6.0 10.0 12.0 0.63 2.77 
18 6.0 10.0 12.0 0.35 3.89 
19 6.0 10.0 12.0 0.83 3.19 
20 6.0 10.0 12.0 0.59 3.17 
21 6.0 10.0 12.0 0.62 2.08 
22 6.0 10.0 12.0 0.80 2.84 
23 6.0 10.0 12.0 0.43 3.08 


2.2 Experiment plots 


The plot designed by Liu et al. (2014a) was selected to obtain different row grades and field slopes 
by simultaneously adjusting screws (Fig. 2). To create seepage conditions, two pipes fixed through 
the plot bottom was used to discharge excess water and then to keep the water level approximately 
1 cm lower than the lowest point of the two conjunction ridges when water was supplied into the 
furrow by pipes banded with gauze. The discharged excess water in the furrow was collected 
through flexible pipes. Seepage and runoff from upper side-slope of the ridge were collected from 
the outlet. 

Brown soil, generated from granite with a high sand content, was used in this experiment. This 
type of soil could easily result in seepage. Table 3 lists the textural information of the soil collected 
from the plow layer. To eliminate the influence of large gravel and plant roots on the experiment 
results, the soil was passed through a 10.0-mm sieve after air dried. The plot was first packed with 
soil at a bulk density of 1.6 g/cm? in four 5-cm layers as the plow sole layers, and then a ridge was 
formed at a bulk density of 1.2 g/cm? as the plow layer. 


Table 3 Textural characteristics of the soil used to pack the experimental plot* 


Gravel Sand Silt Clay Bulk density of plow layer Bulk density of plow sole 
(%) (%) (%) (%) (g/cm*) (g/cm’) 
22,2 71.2 28.1 0.7 0.2 1.6 


Note: *, soil texture was classified according to the standard of U.S. Department of Agriculture (USDA). 
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Fig. 2 Rainfall simulation plot. a, screws for adjusting row grade; b, screws for adjusting field slope; c, pipes 
fixed through the plot bottom for keeping the water level approximately 1 cm lower than the lowest point of the 
two conjunction ridges; d, the lowest point of the two conjunction ridges; e, pipes banded with gauze for providing 
a constant discharge; f, flexible pipes connected to pipes c for collecting excess water from furrow; and g, outlet 
for collecting seepage and runoff from upper side-slope of the ridge. 


2.3 Experiment processes 


The experiments were conducted in the Shandong Provincial Key Laboratory of Soil Conservation 
and Environmental Protection, Shandong, China. The seepage regime was created by supplying 
water (3 L/min) to the furrows, wherein the water level was controlled. In the preliminary 
experiments, seepage discharge from the row sideslopes showed an initial increasing trend and then 
remained at a nearly steady state (Liu et al., 2016). To diminish the effect of seepage discharge 
change on soil erosion, water was supplied for 60 min to create stable seepage conditions. The 
seepage discharge of the final 2 min was collected as a stable seepage flow and used for analysis. 

At the end of the seepage process, the water supply was discontinued, and rainfall simulation 
was performed for 30 min. The rainfall intensity was 39(+0.4) mm/h, and the homogeneity 
coefficient was greater than 0.89 under a trough rainfall simulator fixed with Veejet 80100 nozzles 
(Liu et al., 2014b). The runoff mixed with sediment sample was collected every minute. 

The collected seepage and runoff with sediment samples were weighed immediately. The runoff 
mixed with sediment samples passed through forced-air ovens for 12 h under 105°C. The dried 
sediment yield was then weighed, and the total sediment of a rainfall event was obtained by adding 
the 30 individual samples. 


2.4 Data treatment 


The DPS (data processing system) software package (Tang, 2010) was used to build the second- 
order polynomial regression models and estimate the significance of the regression coefficients. 
Three types of regression models to predict the sediment yield were constructed as follows: 
e Aregression model with the factors of row grade, field slope and ridge heigh. 
e Animproved regression model with the measured seepage discharge and the three factors 
of row grade, field slope and ridge height. 
e Animproved regression model with the predicted seepage discharge and the three factors 
of row grade, field slope and ridge height. 


2.5 Predict performance 


The coefficient of determination (R?) and root-mean square error (RMSE), which represents the 
goodness-of-fit and absolute error measurement, respectively (Legates and McCabe, 1999; Shi et 
al., 2009), were used to evaluate the model performance. Higher R? and lower RMSE values 
indicate a higher prediction accuracy. The equations for R? and RMSE are as follows: 
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2 (Siinsed = S i¢predicted) y 
R? =1 i=1 ; (1) 


n 

2 
> (S icmeasured) ~ Si(mean) ) 
i=1 


= (S icmeasurea) T Sisian) 
RMSE =4|-= , (2) 
n 


where n is the total number of data points and i is the i" number of data point; S is the sediment 
yield (kg); and Smeasured, Spredicted ANd Smean are the measured sediment yield (kg), predicted sediment 
yield (kg) and mean sediment yield (kg), respectively. 


3 Results 


3.1 Seepage model and its influencing factors 


Using the measured seepage discharge (Ysp, L/min) as a dependent variable and the factors (row 
grade (x1 (°)), field slope (x2 (°)) and ridge height (x3 (cm))) as independent variables (Table 2), the 
second-order polynomial regression model was established (Eq. 3). Table 4 lists the results of the 
regression coefficient significance test. Figure 3 illustrates the curve and plot of the measured and 
predicted seepage discharge in the 23 treatments. 


Y,,, = 1.837 — 0.299x, + 0.043x, — 0.169x, + 0.023x,? + 0.13x, 


(3) 
+ 0.002x,x, — 0.002x,x, —0.004x,x, (R? =0.759, P = 0.007) 
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Fig. 3 Curve (a) and plot (b) of the measured and predicted seepage discharge 


Table 4 Results of regression coefficient significance tests for seepage discharge and sediment yield 


Seepage discharge Sediment yield 
Factor 
RC SRC t-Test P RC SRC t-Test P 
RG —0.299 2.31 1.92 0.075 1.80 2.52 2.03 0.062 
FS 0.043 0.41 0.33 0.748 1.55 2.71 2.08 0.057 
RH —0.169 -1.31 0.92 0.375 2.80 3.91 2.65 0.019 
RGxRG 0.023 2.18 3.46 0.004 —0.09 -1.52 2.35 0.034 
FSXxFS 0.000 0.00 0.00 0.998 —0.09 3.24 3.78 0.002 
RHxRH 0.013 2.39 1.93 0.074 0.13 —4.26 3.34 0.005 
RGxFS 0.002 0.19 0.26 0.799 —0.02 —0.29 0.38 0.707 
RGxRH —0.002 0.18 0.19 0.853 —0.03 0.61 0.61 0.549 
FSxRH —0.004 —0.59 0.57 0.581 0.04 1.13 1.05 0.311 
RMSE 0.12 0.67 


Note: RC, regression coefficient; SRC, standard regression coefficient; RG, row grade; FS, field slope; RH, ridge height; x signify the 
interactions between factors; RMSE, root mean square error; —, negative effect of a factor. 
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Eqution 3 explained 75.9% of the variation for seepage discharge at P<0.010. The experiment 
random error accounted for 24.1% of the variation and had a low RMSE value of 0.12. Therefore, 
Equation 3 could be used to predict the seepage discharge. Table 4 showed that the quadratic term 
of the row grade had a significant positive effect on the seepage discharge at P<0.010. The row 
grade had a negative effect, whereas the quadratic term of the ridge height had a positive effect. 
The interactions between the row grade, field slope and ridge height had no significant influence. 
Figure 3 illustrated the curve of the predicted seepage discharge, which fit the observed data well 
(Fig. 3a) and the predicted dots scattered around the observed line (Fig. 3b), although the extent of 
certain values could not be captured precisely. 

Due to the mutual independence of the items in Equation 3, the insignificant items (P>0.100) 
could be ignored as recommended by Tang (2010). Therefore, the simplified Eqution 4 was 
produced, which indicated that the row grade and ridge height might be used to estimate the seepage 
discharge (Ysp) under sufficient water supply conditions and the field slope factor might be omitted. 


Y,» =1.837 — 0.299x, + 0.023x,” + 0.13x,’. (4) 


3.2 Sediment yield model and its influencing factors 


Using a method similar to the seepage discharge model, the regression model for sediment yield 
(Ysy, kg) was obtained (Eq. 5) and tested (Table 4). Equation 5 had a high efficiency coefficient 
(R2=0.743) and low RMSE (0.67), indicating that the sediment yield under seepage conditions 
could be estimated by this equation. Table 4 illustrated the negative effect of the quadratic terms of 
the field slope and ridge height on the sediment yield (P<0.010) as well as the significant negative 
effect of the quadratic term of the row grade (P<0.050). The row grade, field slope and ridge height 
had a positive effect on sediment yield, whereas their quadratic terms had a negative effect, 
indicating that the monofactor effect of these three factors had a converse effect. The interaction of 
the factors had no significant effect on sediment yield even at P<0.100. Based on the P values of 
the items, Equation 5 could be simplified to Equation 6, in which factors and their quadratic terms 
were included and the interactions were eliminated. Equation 6 indicated that the monofactor effect 
could be considered independently by fixing other two factors because no interaction items were 


included in the equation. 
Y,, =—28.38 +1.80x, +1.55x, + 2.80x, —0.09x,’ — 0.09x,” — 0.13x,” 5) 
—0.02x,x, — 0.03x,x, + 0.04x,x, (R? = 0.743, P =0.010) ° 


Yy =—28.38+1.80x, +1.55x, + 2.80x, — 0.09x,” —0.09x,? — 0.13x,’. (6) 


3.3 Improved sediment yield model and its influencing factors 


Equation 7 is the improved regression model for measured sediment yield (Ysy_m; kg) that includes 
the measured seepage (Xm (L/min)) as an independent variable together with the row grade, field 
slope and ridge height. The significance test result of this model is listed in Table 5. The determined 
coefficient for this model was as high as 0.915 with a low P=0.007, meaning a credible estimation. 
The quadratic term of the field slope and interactions between the field slope and seepage discharge 
had a significant negative effect on the sediment yield. Except for the ridge height and seepage 
discharge, all the interactions had a significant effect at P<0.050. When the items with P>0.100 
were eliminated, a simpler form was obtained (Eq. 8). The interaction items were all retained in 
Equation 8, indicating that when the measured seepage discharge was incorporated, all of the 
factors were considered in the sediment yield model. 


Yoy m = 3.87 + 2.69x, +1.25x, — 7.45x, + 78.82x,, — 0.03x," — 0.15x,” 


+0.44x,’ +15.62x ° —0.12x,x, —0.35x,x, + 4.26x,x,, + 0.42x,x, , (7) 
—3.61x,x,, — 7.26x;X,, (R? =0.915, P =0.007) 
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: _= 3.87 + 2.69x, + 78.82x,, — 0,15x,° +15.62x,,” — 0.12x,x, 
— 0.35x,x, + 4.26x,x,, + 0.42x,x, —3.61x,x,, —7.26x,x,, 


Using the predicted seepage discharge (Ysp) from Equation 3 as an input factor, the sediment 
yield (Ysy_p) regression model was built (Eq. 9) with an effective coefficient of 0.893 (P=0.016), 
which suggested a high accuracy of estimation. Table 5 listed the significance test results of the 
regression items. The interaction of row grade and predicted seepage discharge had a positive 
significance effect (P=0.008), which was followed by the quadratic term of seepage discharge at 
P=0.08. Except for these two items, the remaining factors did not have a significant effect even at 
P=0.100. Therefore, the equation might be written as Equation 10, where only the predicted 
seepage discharge and row grade were included. Table 5 and Equations 7 and 9 showed that the 
sediment yield could be more accurately estimated with the measured seepage discharge combined 
than with the predicted seepage discharge incorporated with a lower RMSE (0.38 vs. 0.43) and 
higher R? (0.915 vs. 0.893). 


Ysy p = 1763.99 — 286.24x, + 42.20x, —164.76x, — 943.41x, + 22.14%,’ 
— 0.08x,” +12.71x,” + 38.31x,” +1.78x,x, — 2.13x,x, + 4.62x,X, , (9) 
— 3.862x,x, —1.12x,x, —8.03x,x,  (R° =0.893, P =0.016) 


(8) 


Yoy p =1763.99 + 38.31x,” + 4.62x,x,,, (10) 


where Ysy_p is the sediment yield (kg) with predicted seepage (xp (L/min)). 


Table 5 Sediment yield regression coefficient significance tests incorporated with seepage 


Regression combined with measured seepage Regression combined with predicted seepage 

— RC SRC t-Test P RC SRC t-Test P 

RG 2.69 3.76 2.66 0.026 —286.24 400.22 0.00 1.000 

FS 1.25 2.19 1.72 0.119 42.20 73.76 0.00 1.000 

RH -7.45 —10.41 1.56 0.153 —164.76 -230.37 0.00 1.000 

SR 78.82 14.26 2.19 0.056 -943.41 —148.65 0.00 1.000 
RGxRG —0.03 —0.58 0.75 0.474 22.14 379.86 0.00 1.000 
FSxFS —0.15 5.29 5.80 0.000 —0.08 —2.88 0.00 1.000 
RHxRH 0.44 14.89 1.67 0.129 12.71 429.03 0.00 1.000 
SRxSR 15.62 4.55 2.20 0.056 38.31 9.87 1.97 0.080 
RGxFS —0.12 -2.22 2.77 0.022 1.78 31.59 0.00 1.000 
RGxRH —0.35 —6.72 3.12 0.012 -2.13 —40.34 0.00 1.000 
RGxSR 4.26 6.05 3.05 0.014 4.62 6.22 3.37 0.008 
FSxRH 0.42 10.75 3.19 0.011 -3.86 -98.06 0.00 1.000 
FSxSR -3.61 -8.19 3.37 0.008 -1.12 -2.29 1.00 0.346 
RHxSR -7.26 -22.60 2.23 0.053 -8.03 -23.44 0.97 0.359 
RMSE 0.38 0.43 


Note: RC, regression coefficient; SRC, standard regression coefficient. 


4 Discussion 


4.1 Performance of the improved model 


In previous researches or models (RUSLE and LISEM), the row grade and ridge height were 
considered as subfactors of the contouring coefficient in a contour ridge system (Hessel et al., 2003; 
USDA-ARS, 2008a) or were used to calculate water storage (e.g., WEPP) (Flanagan and 
Livingston, 1995). Because of insufficient data, no direct models had been established between 
these factors and soil loss. Based on the orthogonal rotatable central composite design, five 
gradients of row grade, ridge height and field slope were determined and the regression model for 
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seepage discharge and sediment yield was obtained. Combined with the measured seepage 
discharge, the sediment yield estimation accuracy was clearly improved. The determined 
coefficient (R°) increased from 0.743 to 0.915 and the RMSE decreased from 0.67 to 0.38 for the 
measured seepage discharge combined regression model. In the predicted seepage discharge 
combined model, the R? was increased to 0.893 and the RMSE decreased to 0.43. The result 
indicated that if seepage discharge is available, Equation 7 or its simplified form Equation 8, might 
be used for a better estimation of sediment yield; otherwise, the seepage discharge could be 
predicted by Equations 3 or 4, and then could be used as an input factor for Equations 9 or 10 to 
increase modeling accuracy. Requiring fewer factors, the simplified forms (i.e., Eqs. 4, 8 and 10) 
can be calculated quickly and used more widely when insufficient data are available. 

To assess the improved performance of the combined seepage model, the predicted sediment 
yield and observed values based on the same data set of the regression models were plotted in 
Figure 4. Compared with the regression model that only included the row grade, field slope and 
ridge height (Fig. 4a), the improved model could better predict the extreme value (e.g., treatments 
9 and 10), while the points of the unimproved model were more scattered (Fig. 4b). The improved 
model with measured seepage discharge accurately estimated sediment yield in the latter treatments 
(treatments17—23), i.e., the zero point of the rotatable central composite of the experiment design 
where the predicted value by the seepage incorporated model presented a straight line. The 
predicted error was caused primarily by the repetition of the zero point treatments where the 
seepage discharge and sediment yield should not fluctuate dramatically. However, the improved 
model with the estimated seepage discharge produced better predictions than the model with the 
measured seepage discharge, which was indicated by more compacted points around the predicted 
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Fig. 4 Curve and plot of the measured and predicted sediment yield (a1 and a2 for the regression model; b1 and 
b2 for the regression model combined with measured seepage discharge; c1 and c2 for the regression model 
combined with predicted seepage discharge) 
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line in Figure 4c.This result illustrated that the improved model with the estimated seepage was 
more effective for the treatments around the zero point (treatments 1—16) and that the variance was 
mainly caused by the factor gradients. The difference between the measured and predicted sediment 
yield in the seepage combined model in these treatments might partly be the result of experimental 
error (e.g., soil homogeneity or filling). Therefore, it can be deduced that the improved model with 
predict seepage discharge is able to produce accurate sediment yield estimations. In addition, the 
combination of the simplified Equation 4 for seepage discharge prediction and Equation 10 for 
sediment yield prediction had an advantage that only two factors (row grade and ridge height) were 
required, rendering this approach relatively easy to use. The simple experimental model generally 
provides great assistance for soil conservation and management (Yan et al., 2013). 


4.2 Effects of influential factors on sediment yield 


The influencing factors for soil erosion on sloping land have attracted significant attention (Lal, 
1990; Knapen et al., 2007; Yan et al., 2008; Shi et al., 2013). In contour ridge systems, ridge height 
and row grade together with field slope were considered the main factors influencing soil erosion 
(Wiyo et al., 1999; USDA-ARS, 2008b; Gebreegziabher et al., 2009). The interactions between 
these factors were recently quantified under drainage conditions (Liu et al., 2014a, b). Under the 
seepage conditions, the effects of these factors were deeply investigated in this study. It is 
noteworthy that when the seepage discharge was combined with the regression model, the effect of 
these factors changed significantly (Tables 4 and 5). Table 4 illustrates that sediment yield was 
mainly affected by these three factors and their quadratic terms, while the interactions between 
them had no significant effect. When the measured seepage discharge included, most of the 
interactions increased to a significant level at P<0.050. However, in the predicted seepage 
discharge combined regression model, the significant items were reduced to only the quadratic term 
of seepage discharge and the interaction between row grade and seepage discharge (Table 5). This 
result indicated that under seepage conditions, when considering the inherent relationship between 
seepage discharge and the factors of row grade, ridge height and field slope (i.e., incorporating the 
quadratic polynomial regression model of seepage discharge), the sediment yield prediction model 
could be simplified (Eq. 10). When combining the measured seepage discharge, which was the 
direct result of the row grade, ridge height and field slope and random factors (e.g., measuring error 
and ridge formation error), more interactions between these factors might be detected and more 
variance of the sediment yield may be explained, with the RMSE decreasing and prediction 
accuracy increasing accordingly (Table 5). 

In the measured seepage incorporated regression model, the significant influencing items 
(P<0.050) focused primarily on two categories: the row grade and its interactions with field slope, 
ridge height and seepage discharge; and the quadratic item of field slope and its interactions with 
ridge height and seepage discharge (Table 5). This indicated that under seepage conditions, the 
factors of row grade, field slope and ridge height had primary effect on soil erosion as under 
drainage conditions in RUSLE model (USDA-ARS, 2008b). But the results were different from 
Liu et al. (2014a), who suggested row grade exhibited no significant effect on sediment yield 
because of contour failure. Figure 5 showed the interactions between the paired factors with 
significant effects (P<0.050) on sediment yield when other factors at the central rotated point with 
a zero code level. The interaction between row grade and field slope could be interpreted as having 
a negative effect on sediment yield, with the maximum sediment yield occurring at the highest row 
grade (10°) and the lowest field slope (5°) (Fig. 5a). However, the interaction between row grade 
and field slope exerted no significant effect on interrill and rill erosion under drainage conditions 
(Liu et al., 2014a, b). This indicated that seepage could make the soil erosion process belong to the 
pattern after the critical erosion slope, at which the peak erosion occurs (Jin, 1995; HU and Jin, 
1999; USDA-ARS, 2008b). The decrease of soil erosion as field slope increased may attribute to 
the decrease of detachment power, i.e., the gravity component perpendicular to the slope, when the 
saturated soil had reached its lowest potential resistance under almost the same seepage discharge 
conditions (Liu et al., 2016). Under the highest grade, a ridge height of approximately 12 cm would 
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Fig.5 Significant interactions on sediment yield 


have a lower sediment yield, meaning a better soil conservation (Fig. 5b), which coincided with the 
monofactor analysis result of Liu et al. (2015). A higher seepage discharge would lead to more 
sediment yield in rows with higher grades (Fig. 5c), the reason may be that a higher seepage 
discharge leads to a higher soil erodibility and lower critical shear stress (Nouwakpo et al., 2010) 
and high grade may reduce the soil slope stability along the eroded rills combined with the effect 
of seepage. Field slope and ridge height had a positive interaction on sediment yield, suggesting 
that an increased sediment yield might occur under higher field slopes and ridge height (Fig. 5d), 
which also was suggested by Liu et al. (2014a) under drainage conditions. This might be explained 
as that field slope and ridge height could enhance the soil accumulation slopes leading to more soil 
matrix collapse that had been observed during the erosion process. A higher field slope means a 
lower runoff detachment power, and an increased seepage could lead to a lower erosion resistance 
(Nouwakpo et al., 2010). This two opposite influence might result that the interaction between field 
slope and seepage discharge was negative (Fig. 5e). With a high absolute value of the standard 
regression coefficient (22.6 in Table 5), the interaction between ridge height and seepage discharge 
was important, although the P value was only 0.053. The brightly colored areas in Figure 5f 
indicated a complicated interaction between ridge height and seepage discharge on sediment yield. 
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However, the reason could not be clarified based on the present information. Further study may 
need to reveal the interaction of seepage discharge and ridge height. Base on Figure 5, when one 
factor is fixed, the other factor may be manipulated to produce a minimum sediment yield, and 
subsequent agricultural treatments may be adopted to prevent soil loss. 


4.3 Problems and further studies 


The aggravating effect of seepage on soil erosion was presented by laboratory and field experiments 
and observations (Huang, 1998; Fox and Wilson, 2010; Nouwakpo and Huang, 2012). During the 
water supply period, seepage showed an increasing trend and then remained at an almost steady 
state when the soil was nearly water saturated. Different amounts of seepage discharge could lead 
to different amounts of soil loss (Fox et al., 2007; Nouwakpo et al., 2010). To diminish the effect 
of seepage discharge change on sediment yield, the seepage discharge used as the primary seepage 
condition in this study was measured in the last 2 min of the water supply period. During the rainfall 
simulation period, the seepage discharge may decrease due to void clogging by the detached soil 
particles (Fox and Wilson, 2010), or soil compaction imposed by rainfall and the water infiltration 
(Fohrer et al., 1999), which may affect the sediment yield. However, the seepage discharge mixed 
with the overflow was difficult to be measured separately. Therefore, assessing the effect of seepage 
and its interactions with microtopography (e.g., ridge height) on soil loss during the rainfall event 
remains a difficult task for the further investigation. Further study may need to reveal the interaction 
of seepage discharge and ridge height. 

The hydraulic gradient had an important influence on the critical shear stress, rill erodibility and 
sediment yield (Zheng et al., 2000; Nouwakpo et al., 2010). In a contour ridge system, different 
hydraulic gradients may be created by the water level in the furrows. Although the water level was 
maintained at approximately 1 cm lower than the lowest area of the ridge, the water level could not 
be used as a factor in this study, because it was influenced by field slope and ridge height 
simultaneously. In addition, the water in the furrow was supplied by pipes, but not the runoff 
generated from rainfall. Therefore, determining the effect of rainwater level on the seepage 
discharge and the factors of soil erosion (e.g., critical shear stress and soil erodibility) maybe of 
certain use in revealing the soil erosion process in contour ridge systems. The rainfall intensity and 
regimes play an important role in the loss of soil both individually and interactively with other 
factors (Huang and Laften, 1996; Fang et al., 2012; Liu et al., 2014a). A high rainfall intensity 
might lead to Hortonian overland flow (Ziegler et al., 2007), whereas a low rainfall intensity over 
a long period of time might lead the soil to be saturated with water (Liu et al., 2011), easily resulting 
in seepage. Experiments with different gradients of rainfall intensity also should be addressed in 
soil erosion studies within contour ridge systems under seepage conditions. 


5 Conclusions 


In contour ridge systems, seepage could affect soil erosion. In this study using measured or 
predicted seepage discharge as an input factor, sediment yield estimation models were created to 
increase the R? value from 0.743 to 0.915 and 0.893. Moreover, the improved prediction model 
with the measured seepage discharge was able to detect more significant interactions between the 
influencing factors. When combined with the predicted seepage discharge, the prediction model for 
sediment yield could be substantially simplified. 

Although certain boundary conditions were included, such as the seepage discharge was assumed 
in a nearly steady state and only one rainfall intensity was performed, the results showed that 
considering seepage discharge as an input factor could greatly improve the accuracy of sediment 
yield estimations in contour ridge systems. The effect of rainfall intensities, hydraulic gradients and 
seepage discharge change processes should be addressed in further studies. 
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